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We show that increased plasma superoxide dismutase 1 
(SOD1) levels are statistically significant predictors of the 
failure of pentavalent antimony treatment for cutaneous 
leishmaniasis caused by Leishmania braziliensis. In Leish- 
mania amazonensis-iniected patients, host SOD1 levels can 
be used to discriminate between localized and drug-resistant 
diffuse cutaneous leishmaniasis. Using in situ transcriptom- 
ics (nCounter), we demonstrate a significant positive correla- 
tion between host SOD1 and interferon a/p messenger RNA 
(mRNA) levels, as well as interkingdom correlation between 
host SOD1 and parasite SOD2/4 mRNA levels. In human 
macrophages, in vitro treatment with SOD1 increases the 
parasite burden and induces a diffuse cutaneous leishmania- 
sis-like morphology. Thus, SOD1 is a clinically relevant bio- 
marker and a therapeutic target in both localized and diffuse 
cutaneous leishmaniasis. 
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Leishmaniasis is endemic in several parts of the world, with a 
global prevalence of >12 million cases. Mainly classified in 
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2 distinct clinical forms, leishmaniasis can affect the skin (cuta- 
neous leishmaniasis [CL]) or viscera (visceral leishmaniasis 
[VL]). CL, with 1.5-2 million new cases per year, is an emerging 
infectious disease in several countries, owing to behavioral and 
environmental changes and human immunodeficiency virus 
coinfection [1]. The clinical spectrum of CL includes localized 
CL (LCL), a mild form; diffuse CL (DCL), a disfiguring form; 
and mucocutaneous leishmaniasis (MCL). Leishmania brazil- 
iensis causes LCL and MCL, whereas Leishmania amazonensis 
causes LCL and, sporadically, DCL [2-4]. A low parasite load 
and a variable tendency for self-healing characterize LCL, but 
treatment is strongly recommended, to reduce scarring and pre- 
vent progression to severe MCL. Pentavalent antimonials (Sb v ) 
are currently the first-line treatment for LCL, but because 24.4% 
of the treated patients do not achieve cure with a single cycle of 
Sb v , >2 cycles are indicated to heal the lesions [5]. DCL is char- 
acterized by multiple nonulcerative nodules with an extremely 
high parasite burden, as well as a nonhealing or relapsing phe- 
notype [1,3,4, 5]. DCL is typically refractory to Sb v and second- 
line treatment. Recently, miltefosine treatment induced clinical 
improvement and reduced parasite burden in a cohort of pa- 
tients with DCL, but treatment suspension led to relapse [6]. 

As obligatory intramacrophage parasites, Leishmania species 
have developed several evasion mechanisms to avoid human 
macrophage leishmanicidal activity. One such mechanism 
may involve superoxide dismutase (SOD), a conserved isoen- 
zyme responsible for detoxifying reactive superoxide both in 
Leishmania and human organisms. Previous studies have dem- 
onstrated that SOD-deficient New World Leishmania species, as 
well as Old World Leishmania species, have enhanced sensitiv- 
ity to exogenous superoxide in axenic culture and reduced sur- 
vival in mouse and human macrophages [7, 8]. We have 
previously shown an interferon Pj (IFN-|3) -induced SOD1- 
mediated inhibition of superoxide-dependent parasite killing 
in Leishmania-miected human macrophages in vitro [9]. In ad- 
dition, patients with CL have increased SOD activity in plasma, 
compared with healthy subjects [10]. We hypothesized that 
plasma levels of host SOD1 might reflect clinical status during 
CL and might therefore represent a candidate biomarker. 

PATIENTS AND METHODS 

This study was approved by the Ethics Committee of the Gon- 
calo Moniz Research Center. Informed consent was obtained 
from all patients and healthy controls. CL was diagnosed in pa- 
tients as described elsewhere [11], according to characteristic 
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lesion morphology, positive skin test results, seropositivity to 
Leishmania antigen, and/or presence of parasites in the lesion. 
A total of 58 patients with LCL due to L. braziliensis (27 males; 
mean age [±SD], 29.6 ± 2.3 years) were recruited and treated in 
2 outpatient clinics (Jequie and Tiquirica-BA, NE Brazil) cover- 
ing the same rural area. Ten patients with LCL due to L. ama- 
zonensis (6 males; mean age [±SD], 38.8 ± 6 years) were 
recruited and treated at the Professor Edgar Santos University 
Hospital (Salvador-BA, NE Brazil). Plasma samples from pa- 
tients with LCL were collected at diagnosis (before treatment). 
In addition, a paired sample was obtained from 13 patients at 
the time of definite clinical cure (range, 68-417 days). Patients 
with LCL (due to either L. braziliensis or L. amazonensis) re- 
ceived standard intravenous pentavalent antimony (Glucan- 
time; PJiodia, 20 mg/kg/day for 20 days). Cure was defined by 
the complete scarring of lesions, without induration and relapse 
during 2 years of follow-up. Although all L. amazonensis- 
infected patients were cured within 90 days following a single 
treatment cycle, L. braziliensis-infected patients needed a 
mean (±SD) of 2.2 ± 0.2 treatment cycles over 190.5 ± 16.0 
days to achieve cure, with treatment for 23 of 58 patients 
(39.7%) classified as failing (>2 cycles). Eight patients with 
DCL (3 males; mean age [±SD], 31.5 ±7.1 years) infected 
with L. amazonensis were recruited and treated at the Presidente 
Dutra University Hospital-HUPD (Sao Lufs-MA, NE Brazil). 
Patients with DCL had a mean disease duration (±SD) of 
10.9 ±2.1 years, during which several different therapeutic 
schemes (mean number [±SD] per patient, 2.7 ± 0.9, including 
liposomal amphotericin B and combinations of Glucantime, in- 
terferon y, and pentamidine) were used, with no or only a tran- 
sient clinical effect. Since DCL is a rare clinical manifestation, 
no samples from untreated patients (at diagnosis) were available 
for comparison. SOD1 plasma levels were quantified using a 
human SOD1 enzyme-linked immunosorbent assay kit (Cal- 
biochem). RNA was extracted from lesion biopsy specimens 
from patients with LCL and those with DCL, as well as skin bi- 
opsy specimens from healthy controls, using Trizol, followed by 
an additional purification using RNeasy (QIAgen Benelux, 
Venlo, the Netherlands). Using in situ transcriptomics, we si- 
multaneously quantified host and parasite RNAs in lesion biop- 
sy specimens from patients with LCL and those with DCL, as 
well as skin biopsy specimens from healthy controls, by 
nCounter technology (NanoString, Seattle, WA), based on mo- 
lecular bar coding of target RNA transcripts and digital detec- 
tion at the femtomolar range, using direct hybridization, 
without reverse transcription or amplification [12]. Human 
SOD1, SOD2, SOD3, IFNA1, IFNA2, IFNA4, IFNB1, and 
L. braziliensis and L. amazonensis SOD 1-5 messenger RNA 
(mRNA) levels were quantified in situ, in addition to mRNAs 
of several housekeeping genes (GUSB, G6PD, GAPDH, and 
HPRT1), for normalization, as well as leukocyte-specific genes 
(CD3, CD14, CD19, CD56, and CD45). Quantification of the 



parasite burden in vitro was performed as described elsewhere 
[9]. Briefly, monocyte-derived human macrophages were infect- 
ed with L. amazonensis (MHOM/BR/87/BA125) and treated 
with SOD1 (Sigma) for 48 hours, followed by extensive washing, 
staining with hematoxylin/eosin, and counting of intracellular 
amastigotes (100 cells, with duplicate analysis for each sample). 
Parametric and nonparametric tests were performed according 
to Kolmogorov-Smirnov test for normality. For multiple com- 
parisons, the Kruskal-Wallis test with the Dunn posttest and 1- 
way analysis of variance with a posttest for linear trend were 
used. For comparison between 2 groups, the F test, the 
Mann-Whitney U test, the t test, or the Wilcoxon test was 
used. To identify biomarkers, Receiver Operating Characteristic 
(ROC) curve analysis was used. All tests were 2 tailed, and dif- 
ferences were considered significant at P values of <.05. 

RESULTS 

We measured plasma levels of SOD1 in healthy controls, pa- 
tients with LCL (due to L. braziliensis or L. amazonensis), 
and patients with DCL (due to L. amazonensis) and investigated 
its potential role as a biomarker. The mean SOD1 level (±SEM) 
among patients with LCL due to L. braziliensis (90.6 ± 8.0 pg/ 
mL) was 3 -fold greater than that among healthy controls 
(30.2 ± 3.9 pg/mL; P < .01, by the Kruskal-Wallis; P < .001, by 
the Dunn posttest), whereas the mean level (±SEM) among pa- 
tients with LCL due to L. amazonensis (93.1 ± 24.7 pg/mL) was 
3.1-fold greater than that among healthy controls (P< .01, by 
the Kruskal-Wallis test; P< .01, by the Dunn posttest). On the 
other hand, mean SOD1 levels (±SEM) among patients with 
DCL (407.3 ± 88.0 pg/mL) were 4.4-fold greater than those 
among patients with LCL due to L. amazonensis (P = .0005, 
by the Mann-Whitney U test; Figure 1A). To investigate wheth- 
er increased SOD1 levels during LCL is related to active disease, 
the SOD1 level was also measured in paired samples of from 13 
patients with LCL due to L. braziliensis before treatment (at di- 
agnosis) and after treatment (at definite clinical cure). SOD1 
levels decreased 1.6-fold between pretreatment and posttreat- 
ment samples (99.2 ± 20.81 pg/mL vs 60.4 ± 5.2 pg/mL; 
P = .027, by the Wilcoxon matched pairs test; Supplementary 
Figure 1A). Since variances differed significantly between con- 
trols and patients with LCL due to L. braziliensis (P < .0001, by 
the F test), we hypothesized that heterogeneous SOD1 levels in 
patients with LCL might reflect differences in clinical status 
and/or therapeutic response. SOD1 levels did not significantly 
correlate with age, sex, disease duration, lesion size, or healing 
time among patients with LCL (data not shown). However, 
SOD1 plasma levels were significantly associated with therapeu- 
tic response, quantified as the number of treatment cycles re- 
quired for cure (P= .0069, by the Kruskal-Wallis test; P< .001, 
by the Dunn posttest; Figure 1C). By use of ROC curve analysis, 
successful (1 cycle) and failing (>2 cycles) therapeutic response 
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Figure 1. Superoxide dismutase 1 (S0D1) is a biomarker for therapeutic failure in cutaneous leishmaniasis. A, S0D1 measurements in plasma samples 
from 20 healthy controls, 58 patients with localized cutaneous leishmaniasis (LCL) due to Leishmania braziliensis), 10 patients with LCL due to Leishmania 
amazonensis, and 8 patients with diffuse cutaneous leishmaniasis (DCL) due to L amazonensis. Box plots represent the median and interquartile range of 
each group. P< .001, by the Kruskal— Wallis test; and ***p< .001, by the Dunn posttest, for healthy controls vs patients with LCL due to L braziliensis; 
**P< .01, by the Dunn posttest, for healthy controls vs patients with LCL due to L. amazonensis; and ***p< .001, by the Mann-Whitney fy test, for patients 
with LCL due to L. amazonensis vs patients with DCL. B, Plasma S0D1 level significantly discriminates patients with LCL due to L. amazonensis from 
patients with DCL due to L. amazonensis (area under the receiving operating characteristic [ROC] curve, 0.95; ***P=.001). C, S0D1 measurements in 
plasma samples from patients with LCL due to L. braziliensis, classified according to number of treatment cycles (P= .0069, by the Kruskal— Wallis test; 
**P< .001 , by the Dunn posttest). Three outliers (by the Tukey test) are shown as single dots. D, Plasma S0D1 level significantly predicts therapeutic 
response in patients with LCL due to L. braziliensis (success is defined as achievement of cure after 1 treatment cycle; failure is defined as lack of 
cure after >2 treatment cycles; area under the receiving operating characteristic curve, 0.83; ***p= .00025, without outliers, as in panel C). 



could be significantly discriminated through SOD1 plasma lev- 
els (area under the ROC curve, 0.83; P = .00025), authenticating 
its clinical utility as a biomarker (Figure ID). Above a cutoff of 
122 pg/mL, 9 of 22 patients (40.9%) in whom treatment was 
failing could be identified with 100% specificity (ie, 0 patients 
cured with 1 treatment cycle displayed SOD1 levels of >122 
pg/mL). Strikingly, at a similar cutoff of 112 pg/mL, patients 
with DCL could be discriminated from treatable patients 
with LCL (due to L. amazonensis) with 100% sensitivity and 
90% specificity (area under the ROC curve, 0.95; P=.001; 
Figure 15). 



To test whether these high plasma SOD1 levels in patients 
with CL might be directly related to an increased parasite bur- 
den, human macrophages from healthy donors were infected 
with L. amazonensis (5:1) and cultivated in the absence or pres- 
ence of recombinant SOD1 protein. As shown in Figure 2A, 
SOD1 treatment significantly increased parasite load (P= .042, 
by the t test). Moreover, infected human macrophages treated 
with recombinant SOD1 protein presented cytomorphological 
features (Figure 2B) strikingly similar to the vacuolized and 
highly parasitized macrophages from cutaneous lesions of pa- 
tients with DCL (Figure 2C). 
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Figure 2. Superoxide dismutase 1 (S0D1) correlates to parasite burden in vitro and in situ in diffuse cutaneous leishmaniasis (DCL). A and B, Human 
macrophages were infected with Leishmania amazonensis (5:1 ratio) for 4 hours and then treated in the absence or presence of S0D1 protein (1 75 U/mL) for 
48 hours. Cells were fixed on glass slides and stained with hematoxylin and eosin. A, Parasite burden was quantified; each bar represents the mean ± stan- 
dard error of the mean for 3 donors (*P= .042, by the ftest). B, Infected macrophages (1000x magnification) left untreated (left panel) or treated with S0D1 
(right panel). C, Cutaneous lesion biopsy specimen from a representative patient with DCL stained with hematoxylin and eosin (1000x magnification). D, 
Positive correlation between human S0D1 and parasite S0D2 (r = 0.98, P= .019) and S0D4 (r = 0.99, P= .0026) messenger RNA (mRNA) levels (normalized 
to G6PD), quantified by nCounter in situ in biopsy specimens from patients with DCL (n = 4). 



Finally, using in situ transcriptomics (nCounter), we demon- 
strated a significant positive correlation between host SOD1 
mRNA level and parasite SOD2 (P= .019, r = 0.98) and SOD4 
(P = .0026, r = 0.99) mRNA in DCL. Taken together, these data 
suggest that systemic SOD1 might play a causal role in both the 
clinical and anatomopathological phenotype of DCL. 

DISCUSSION 

Up to now, no noninvasive biomarker for therapeutic response 
to Sb v has been described in CL. In the present work, we dem- 
onstrate host SOD1 as a biomarker for therapeutic failure of 
first-line Sb v in New World patients with CL. SOD1 plasma 
levels could significantly discriminate between successful and 
failing therapeutic response in patients with LCL (due to 
L. braziliensis) and between treatable patients with LCL (due 
to L. amazonensis) and refractory patients with DCL. Quantifi- 
cation of SOD1 plasma level at diagnosis represents a rapid, 



sensitive, and inexpensive new tool in the clinical management 
of CL. Further research will be necessary to test its predictive 
value in CL caused by other Leishmania species, since in Peru 
failure of antimonial therapy is significantly associated with in- 
fections due to L. braziliensis and Leishmania peruviana but not 
with infection due to Leishmania guyanensis [5]. By identifying, 
at the time of diagnosis, up to 40.9% of patients (with 100% spe- 
cificity) in whom antimonial therapy will fail (Figure 1C), sec- 
ond-line treatment (eg, amphotericin B or miltefosine) could be 
provided as a first choice, thereby diminishing public health 
costs and avoiding unnecessary and prolonged exposure of 
the patients to the often severe side effects of intravenous pen- 
tavalent antimony therapy. 

Although SOD1 is an intracellular antioxidant enzyme, our 
results suggest a direct role for increased systemic (plasma) 
SOD 1 as a trigger of parasite replication in infected macrophages. 
Elsewhere, treatment of infected human macrophages with pu- 
rified SOD1 in vitro led to strongly increased parasite burden 
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and the appearance of large parasitophorous vacuoles, reminis- 
cent of typical nodular lesions in DCL [13]. Moreover, this cor- 
relation between SOD1 and parasite burden was recapitulated 
in situ in DCL lesions, as evidenced by a strong linear relation- 
ship (r = 098-0.99) between host (human SOD1) and parasite 
(L. amazonensis SOD2 and SOD4) SODs. This parallel cross- 
regulation at the transcriptional level between superoxide- 
quenching enzymes that are hundreds of millions of years 
distant in evolution is a striking example of interkingdom sig- 
naling. To our knowledge, this is the first description of inter- 
kingdom signaling in protozoan-animal interaction. In keeping 
with our cytomorphological data (Figure 25), this phenomenon 
is highly specific for cytosolic SOD1, since cross-kingdom 
correlations were not observed for mitochondrial SOD2 or extra- 
cellular SOD3 with any of the parasite mRNAs. The extremely 
long (mean disease duration [±SD], 10.9 ±2.1 years) uncon- 
trolled parasite replication in patients with DCL might have led 
to selection of SOD-overexpressing parasite clones [7, 8] . Alterna- 
tively, increased systemic levels of SOD1 in patients with LCL and 
those with DCL might reflect a parasite escape mechanism from 
host leishmanicidal activity mediated by type I IFN, as suggested 
by our previous work [9] and recently confirmed in murine and 
human macrophages [14]. In agreement with this hypothesis, we 
found a significant positive correlation between SOD1 mRNA 
and total IFN-a/p mRNA levels, using in situ transcriptomics 
(nCounter; Supplementary Figure 2A). 

Finally, since all of our patients with DCL also experienced 
failure of multiple other treatment schemes (including pentami- 
dine, amphotericin B, and IFN-y, alone or in combination), an 
increased plasma SOD1 level truly reflects broad therapeutic fail- 
ure. The dramatic increase in DCL also reinforces the idea that 
SOD1 is not a mere surrogate marker but a key molecule in 
DCL pathogenesis, driving perennial parasite replication and sys- 
temic dissemination, as suggested by Figure 2A-D. Hence, SOD1 
is also a candidate therapeutic target in CL, reinforcing our pre- 
vious suggestion of pharmacologic SOD1 inhibition by DETC 
[15] as a viable therapeutic alternative in LCL and DCL in par- 
ticular, for which no effective treatment options exist. 

In conclusion, our results indicate that host SOD1 is a clini- 
cally useful biomarker for therapeutic failure, as well as a pos- 
sible therapeutic target, in LCL and DCL. 
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